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ABSTRACT 

In TeV scale gravity models, for dilepton production at hadron colliders, we present 
the NLO-QCD corrections for the double differential cross section in the invariant 
mass and scattering angle. For both ADD and RS models, the quantitative impact 
of QCD corrections for extra dimension searches at LHC and Tevatron are investi- 
gated. We present the K-factors for both ADD and RS models at LHC and Tevatron. 
Inclusion of QCD corrections to NLO stabilises the cross section with respect to scale 
variations. 
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1 Introduction 



Extra dimension scenarios are now essential part of the studies of physics beyond 
the Standard Model (SM). They provide an alternate view of the hierarchy between 
the electroweak and the Planck scale. Some of these extra dimension models in- 
voke the brane world scenarios to hide the extra spacial dimensions from current 
observation. Two such models that are phenomenologically widely studied are the 
Arkani-Hamed, Dimopoulos and Dvali (ADD) [1] and the Randall-Sundrum (RS) 
[2] models. 

In the ADD case the compactified extra dimensions could be large and the large 
volume of the compactified extra spacial dimension would account for the dilution 
of gravity in 4-dimensions and hence the hierarchy. In this model, new physics 
can appear at a mass scale of the order of a TeV. A viable mechanism to hide the 
extra spacial dimension, is to introduce a 3-brane with negligible tension and lo- 
calise the SM particles on it. Only the graviton is allowed to propagate the full 
4 + d dimensional space time. As a consequence of these assumptions, it follows 
from Gauss Law that the effective Planck scale M P in 4-dimension is related to the 
4 + d dimensional fundamental scale M s through the volume of the compactified 
extra dimensions [1]. The extra dimensions are compactified on a torus of common 
circumference R. The number of extra spacial dimension possible is d > 2 from 
current experimental limits on deviation from inverse square law [3]. The space 
time is factorisable and the 4-dimensional spectrum consists of the SM confined to 
4-dimensions and a tower of Kaluza-Klien (KK) modes of the graviton propagating 
the full 4 + d dimensional space time. 

The interaction of the KK modes hffi with the SM fields localised on the 3-brane 
is given by 

1 oo 

C int ^- w ^T^(x)h^(x), (1.1) 
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where T^ v is the energy-momentum tensor of the SM fields on the 3-brane. The zero 
mode corresponds to the usual 4-dimensional massless graviton. The KK modes 
are all M P suppressed but the high multiplicity could lead to observable effects at 
present and future colliders. The Feynman rules are given in [4, 5]. 

In the RS model there is only one extra spacial dimension and the extra dimen- 
sion is compactified to a circle of circumference 2L and further orbifolded by iden- 
tifying points related by y — > —y. Two branes are placed at orbifold fixed points, 
y = with positive tension called the Planck brane and a second brane at y = L with 
negative tension called the TeV brane. For a special choice of parameters, it turns out 
that the 5-dimensional Einstein equations have a warped solution f or < y < L with 
metric g fJiU (x p ,y) = exp(—2ky) r)^, g m = and g yy = 1. This space is not factoris- 
able and has a constant negative curvature — AdS 5 space-time, k is the curvature of 
the AdS 5 space-time and r]^ is the usual 4-dimensional flat Minkowski metric. In 
this model the mass scales vary with y according to the exponential warp factor. If 
gravity originates on the brane at y — 0, TeV scales can be generated on the brane 
at y = L for kL ~ 10. The apparent hierarchy is generated by the exponential warp 
factor and no additional large hierarchies appear. The size of the extra dimension is 
of the order of Mp 1 . Further it has been showed that [6] the value of kL can be sta- 
bilised without fine tuning by minimising the potential for the modulus field which 
describes the relative motion of the two branes. In the RS model graviton and the 
modulus field can propagate the full 5-dimensional space time while the SM is con- 
fined to the TeV brane. The 4-dimensional spectrum contains the KK modes, the 
zero mode is M P suppressed while the excited modes are massive and are only TeV 
suppressed. The mass gap of the KK modes is determined by the difference of the 
successive zeros of the Bessel function J\{x) and the scale m = k e~ wkL . As in the 
ADD case the phenomenology of the RS model concerns the effect of massive KK 
modes of the graviton, though the spectrum of the KK mode is quite different. 
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In the RS model the massive KK modes hffi (x) interacts with the SM fields 

i p TvkL oo 

An* ~ - W T^{x)h^{x) - — ^T^(x)h^(x) , (1.2) 

where T^ v is the energy-momentum tensor of the SM fields on the 3-brane at y = L. 
The masses of hffi (x) are given by M n = x n k e~ nkL , where x n are the zeros of the 
Bessel function Ji(x). In the RS model there are two parameters which are c = 
k/Mp, the effective coupling * and Mi the mass of the first KK mode. Expect for an 
overall warp factor the Feynman rule of RS is the same as those of the ADD model. 

Next to leading order (NLO) QCD corrections have been recently calculated in 
the ADD case for e + e~ — > hadrons [7] and various distributions of invariant lepton 
pair production at both LHC and Tevatron [8]. This was further extended to the RS 
case [9]. In this paper, for the ADD and RS models, we consider the un-integrated 
distribution with respect to cos 9* to NLO in QCD, where 9* is the scattering angle of 
the lepton with an initial hadron in the c.o.m frame of the lepton pair. This is partic- 
ularly important in the dilepton production case to achieve maximum sensitivity to 
the model parameters, as cos 9* integrated cross section is independent of the inter- 
ference between SM and gravity [8]. To leading order (LO), this double differential 
da I dQ I d cos 9* was analysed in [10]. At hadron colliders the NLO-QCD corrections 
are important especially in models of extra dimension as gluon-gluon subprocess 
contributes at the same LO as quark-antiquark subprocess. DO Collaboration re- 
cently reported searches for large extra dimensions in the dimuon channel for the 
double differential cross section [11], this updates the earlier Run-I results [12]. The 
first direct search of the RS KK modes using the dileptons have been reported by 
DO Collaboration [13]. 

Rest of the paper is organised as follows: In section 2 we evaluate the NLO 
coefficient functions to the subprocess that contribute to the double differential cross 
section da/dQ/d cos 9*. Finally in section 3, we discuss the impact of the NLO results. 
*An alternate definition is c = k/M P , where M P = M p /\/&k, hence c = c \f&ir. 
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2 Drell-Yan cos 6* distribution 



We consider the Drell-Yan process and study the double differential cross section 
with respect to the invariant mass of the final lepton pair and cos 9* the cosine of 
the angle between the final state lepton momenta and the initial state hadron in the 
c.o.m frame of the lepton pair § . The relevant kinematical formulation is detailed in 
[8]. In the QCD improved parton model, the hadronic cross section can be expressed 
in terms of partonic cross sections convoluted with appropriate parton distribution 
functions. The coefficient functions to NLO in QCD are evaluated for both ADD and 
RS models. The difference between the two models depend on the spectrum of the 
KK modes and hence summation of the KK modes that contribute to the dilepton 
production leads to different results [8, 9]. 

The hadronic part involves the computation of various processes that contribute 
to Q or X F or rapidity distributions that are presented in the reference [8]. The 
angular distributions which are "odd" in cos 9* come mainly from the interferences 
terms. The non-vanishing odd contribution in the standard model sector comes from 
the interference of photon mediated processes with Z-boson mediated processes. 
We also find that non-vanishing odd contributions come from the interference of 
standard model diagrams with the graviton exchange diagrams. These inference 
diagrams are absent in the computation of Q, X F and rapidity distributions where 
only even functions of cos 9* contribute. We have regularised all the divergences 
using dimensional regularisation. The mass singularities are removed by the mass 
factorisation, for details refer to [8]. 

We first present the angular distribution which is "even" in cos 9*. 

da PlP2 f 1 f 1 f 1 

2S ^T7 ^ = Y,^SM, q dxi / dx 2 / dz 5(t - zx x x 2 ) 

dQ 2 d cos 9* g Jo Jo Jo 

§ An alternate definition of the angle has been considered in [14] to study the lepton helicity dis- 
tribution in polarised Drell-Yan process. 
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X 



H qq {x u x 2 , &) (Af^ /Z (z, Q 2 , &) + a s A q f j/Z (z, Q 2 , &) 



+H qg ( Xl ,x 2 , fi 2 F )a s A ( q f^ z (z, nl) 
+H gq (x u x 2 ,fi 2 F )a s A^ z (z,fi 2 F ) 



+ J2 :F SMGR,q J q dx l J Q dX2 f dzS(T-ZX!X 2 ) 



X 



H qq ( Xl , x 2 , (A q f G " /z (z, Q 2 , &) + a s A q f G " /z (z, Q 2 , & 



+H qg (x 1 ,x 2 ,^ F )a s A q 1 J' G ^ z (z,^ F ) 
+H gq (x 1 ,x 2 ,fx 2 F )a s A g 1 q ^/ z (z, f x 2 F ) 



^ J~gr i dx\ / dx 2 / dz 5(t — zx\X 2 ) 
„ Jo Jo Jo 



X 



H qq (x u x 2 , /4) (A g f G (z, Q 2 , + a s A q f G (z, Q 2 , &)) 



+H qg (x 1 ,x 2 , fi 2 F )a s A q f G (z, Q 2 , /4) 
+H gq ( Xl ,x 2 , fj, 2 F )a s A g ^ G (z, Q 2 , /4) 
+H gg ( Xl ,x 2 , (A g f G (z, Q 2 , /4) + a s A g ^ G (z, Q 2 , ^)) 



(2.1) 



where H ab (xi,x 2 , fi 2 F ) are the renormalised partonic distributions and A ab (z, Q 2 , fj, 2 F ) 
are the coefficient function corresponding to various subprocess cross section to 
NLO in QCD. The factors TsM, q , Pgr correspond to pure SM and gravity (GR) part 
respectively and are given in [8], the factor that corresponds to the interference of 
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SM and gravity is 



JSMGR,q = TI~V- ) \Q )\ 



4:71 



Q\Q 2 -M 2 ) 



(2.2) 



((Q 2 - M\Y + M 2 Y\) clsl 

where a is the fine structure constant, k = y/16n/Mp. The summation of the KK 
modes leads to V(Q 2 ) for ADD and RS case has been given in [8, 9]. 
The leading order results read 
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and the next to leading order results read 
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where C F = {N 2 - 1)/2N, C A = N and T> = 1/2 are the SU{N) colour factors and 
n/ is the number of flavours. The "plus" functions appearing in the above results 
are the distributions which satisfy the following equation 



where 



r dz u 

Jo 



/+(*) 



£dzf(z)(g(z)-g( 



4n'(l - z) 
1-z 



0,1 



+ 



and g(z) is any well behaved function in the region < z < 1. 

In Eqs. (2.3) and (2.4) the term (1 — 3cos 2 #*) corresponds to the interference 
between the SM and GR and within the SM interference between 7 and Z dia- 
grams. Though this combination is even in cos 9* it vanishes in the angular inte- 
grated cross section and also does not contribute to the forward-backward asymme- 
try A FB - Hence the un-integrated cross section is very useful to study this contribu- 
tion to the interference effect in the Drell-Yan process. 
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We present below the angular distribution which is "odd" in cos 9* 

da™ _ 

q* 

1 



dQ 2 d cos 9* 



y2SFsM, q ( dxi ( dx 2 [ dz5(T-ZXiX 2 ) 
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(2.5) 



The constants 5F S M,q, SFsMGR, q are given by 
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(2.6) 
(2.7) 



The renormalised incoming partonic fluxes are defined by 

SH qq ( Xl ,x 2 , //|) = fq'ixufil) f q D2 {x 2 , fi 2 F ) - f^ix^nl) f q 2 {x 2 , fi 2 F ) , 

SHggfaxtoiip) = fg'ixi,^) (fg 2 {x2,H 2 F ) ~ fq 2 {x2^ 2 F )) , 
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6H qg (x 1 ,x 2 ,ii F ) = 8H gq (x2,xi,iJ? F ) . 
The LO coefficient functions corresponding to Eq. (2.5) are 
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(2.10) 



These coefficient functions Eq. (2.9) and (2.10), which are odd in cos#*, are due to 
the interference of 7 and Z in SM and between SM and GR in the full theory. Note 
that the q — > g in the subprocess leads to a negative sign which has been taken 
care of in the renormalised parton fluxes Eq. (2.8). A FB picks up this odd parts 
which contributes to the interference terms. To NLO the A FB coefficient functions 
have been evaluated in [8] 1 and the effects analysed for the ADD case. In the next 
section, the impact of the NLO-QCD correction derived in this section is discussed. 



11 In [8] the last three equations of Eq. (6.11), the RHS should read 
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3 Discussions 



In this section, the effect of the NLO QCD corrections on the angular distribution of 
lepton pair are presented. We present these distributions for the LHC (y/S = 14 TeV) 
and Run II of Tevatron (y/S = 1.96 TeV) for typical values of ADD and RS model 
parameters. The effort here is mainly to emphasise the impact of QCD correction on 
the bounds rather than to extract bounds on Ms. 

For ADD model, we choose the parameters M s = 2 TeV and d = 3. For RS we 
choose c = 0.01, M l = 1.5 TeV(for LHC) and M x = 700 GeV(for Tevatron). The 
SM parameters which enter our analysis are a = 1/137.03604, M z = 91.1876 GeV, 
T z = 2.4952 GeV and sin 2 9 W = 0.227. For the parton density sets, we adopt in 
leading order, the MRST 2001 LO (A = 0.1670 GeV) and in next-to-leading order, the 
MRST 2001 NLO (A = 0.2390 GeV). The renormalisation scale fi R and factorisation 
scale jijp are taken to be equal to Q unless mentioned otherwise. 

For the coefficient functions which are even in cos 9*, the parton density com- 
binations are even under the interchange of x 1 and x 2 , while for the odd terms in 
cos 9*, the parton density combinations are odd under this exchange. Hence, the 
quark-antiquark initiated contributions from these odd terms to LHC cross sections 
are zero, but small contribution comes from quark-gluon initiated processes. This is 
not the case for Tevatron. 

In the SM part, at LO level, the quark-antiquark initiated processes behave as 
1 + cos 2 6* for pure 7 and Z intermediate states, and as cos 6* for 7Z interference 
terms. In the Gravity part, at LO, gluon-gluon initiated process is of the form 1 — 
cos 4 6>*, quark-antiquark process is of the form 1 — 3 cos 2 9* + 4cos 4 #*. The inter- 
ference between SM and GR always behaves as cos 9*, cos 3 9* and 1 — 3 cos 2 9*. At 
NLO, quark-gluon initiated processes contribute to both SM, GR and the interfer- 
ence terms. 

We first discuss the phenomenology at LHC using ADD model. In Fig. la, we 
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plot the angular distribution at Q = 700 GeV with NLO corrected cross sections. We 
find that the gravity contribution is large compared to that of SM. Since the gluon 
flux is large at LHC, the gluon-gluon initiated subprocess dominates over the rest. 
Also, the interference between SM and Gravity is negligible over the entire range of 
cos 9* at this Q = 700 GeV. In the Fig. lb, we have plotted the K-factor at LHC. In 
general K-factor is defined as 



K 



dQd cos 9* 



d<y I NLO (Q,cos9* 
dQd cos 9* 



(3.1) 



where / = SM, TOT, SM means Standard Model, TOT means sum of SM and Grav- 
ity contributions. Since the gluon-gluon initiated process dominates over the rest 
the K-factor is around 1.4 to 1.5 in the entire range of cos#*. In the Fig. lc we have 
plotted the R-ratio in order to check whether the NLO results improve the scale un- 
certainty. Here we have chosen fi R = fi F = fi and we follow the same throughout 
our analysis. The R-ratio is defined as 



rI lo - 



Rnlo 



dQd cos 9* 



d(T NLo(^ = t* o)) 
dQd cos 9* 



dQd cos 9* 



- -1 


d^NLO^)) 




dQd cos 9* 



J Q=700 GeV 



Q=700 GeV 



(3.2) 



(3.3) 



We have chosen 9* = 45° for the plot. As we can see the NLO results improve the 
scale uncertainty. 

Let us now repeat the similar study for the RS model at LHC energies. In the 
Fig. 2a, we have plotted the angular distribution at first resonance M 1 = 1.5 TeV and 
c = 0.01. We find that the gravity contribution is well above the standard model 
one. In particular, the gluon-gluon initiated contribution is the dominant one. Since 
the SM and interference contributions are of the same order and are extremely small 
due to large Q which is 1.5 TeV, the total contribution is purely due to the gluon- 
gluon initiated process. Unlike the ADD case (see Fig. la), the total contribution 
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mainly comes from the gravity mediated process at the first resonance. In Fig. 2b, 
we have plotted the K-factor defined in Eq. (3.1) at the first resonance. Because of 
this, the dominant contribution comes from the gravity mediated processes. Since 
Q — 1.5 TeV, both quark-antiquark as well as the gluon initiated processes contribute 
at the same level because their partonic fluxes are comparable at this energy. The 
shape of the K-factor in Fig. lb looks different from Fig. 2b because at Q = 0.7 
TeV, only gluon initiated process dominates. In Fig. 2c, we have plotted the R-ratio 
defined in Eq. (3.2, 3.3) for RS resonance M 1 = 1.5 TeV, c = 0.01 at 9* = 45°. From 
the plot it is clear that the NLO corrections improve the scale uncertainty. 

Next we discuss the phenomenology at Tevatron. We start with ADD for the 
parameters M s = 2 TeV, d — 3 and Q = 400 GeV In Fig. 3a we have plotted angular 
distribution with NLO results. We find that the SM dominates over the rest. In Teva- 
tron, both the even and odd in cos 0* contribute significantly unlike in the LHC. This 
leads to the asymmetry in the angular distribution as shown in Fig. 3a at Tevatron. 
At Q = 400 GeV the dominant contribution at Tevatron is from SM. In Fig. 3b, we 
have plotted the K-factor at Q = 400 GeV. Since it is the quark-antiquark initiated 
process that dominates, the K-factor is similar to the SM value which is around 1.3. 
Fig. 3c shows the sensitivity of the results to the scale variation. As expected NLO 
improves the result. 

We now study the phenomenology at Tevatron in the case of RS model. We 
choose the RS model parameters Mi = 700 GeV and c = 0.01 which are not ex- 
cluded by the recent searches by D0[13]. In Fig. 4a, we have plotted the angular 
distribution at the first resonance using NLO corrected cross sections. Being in the 
resonance region, gravity mediated process dominates over the SM. At Q = 700 
GeV, the quark initiated processes contribute significantly, as the qq flux dominates. 
In Fig. 4b, we have plotted the K-factor at Q = 700 GeV. Since the quark initiated 
process dominates, the K-factor is close to the SM value. Fig. 4c shows the sensitiv- 
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ity of the results with respect to scale. One can easily notice that NLO gives reliable 
predictions. 

In summary, we have computed the cross sections da / dQ / d cos 9* up to next to 
leading order in QCD. Along with the standard model results, we have presented 
the contributions from all the subprocesses that are due to the graviton in the context 
of TeV-scale gravity models. Our main conclusion is that the NLO QCD corrections 
are very significant at the LHC because of the large incident gluon flux. At the Teva- 
tron where the gluon flux is small, the NLO effects are moderate for ADD and RS in 
the resonance region. But, significantly, at both the colliders the inclusion of the NLO 
QCD corrections help stabilise the cross-section with respect to scale variations. The 
extraction of bounds from both the colliders will, therefore, require the inclusion of 
these NLO QCD corrections. 
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Figure Caption 



Figure 1. (a) The double differential cross section da /dQ/d cos 9* is plotted as a func- 
tion of cos 6* for Q = 700 GeV at LHC. The typical ADD parameters chosen are 
M s = 2 TeV, d — 3. (b) The corresponding K-factor for cos 9* distribution SM and SM 
plus gravity (TOT), (c) Scale variation at LO and NLO as defined in Eq. (3.2), (3.3) 
for Q = 700 GeV and 6* = 45°. 

Figure 2. (a) The double differential cross section da /dQ/d cos 6* is plotted as a func- 
tion of cosO* for Q = 1.5 TeV at LHC. The RS model parameters are M 1 = 1.5 TeV 
and c = 0.01. (b) The K-factor for the distribution in (a) is plotted for the cos0* 
range [-1,1]- (c) The scale variation of the ratio R is plotted as a function of fi/fi at 
the first RS KK resonance region and 9* = 45°. 

Figure 3. (a) For Tevatron energies, da /dQ/d cos 8* is plotted as a function of cos 9* 
at Q = 400 GeV for typical value of ADD parameters M s = 2 TeV and d = 3. (b) The 
K-factor for cos 6* distribution for the same ADD parameters in (a) is plotted, (c) The 
variation of the R-ratio with respect to the scale /i//i for the ADD parameters in (a) 

at 6* = 45°. 

Figure 4. (a) The double differential cross section da /dQ/d cos 6* is plotted as a func- 
tion of cos#* for Q = 700 GeV at the Tevatron. RS parameters M 1 = 700 GeV and 
c = 0.01. (b) The K-factor for the distribution in (a) is plotted for the cos#* range, 
(c) The scale variation of the ratio R is plotted as a function of /i//i for 6* = 45° and 
Q = 700 GeV. 
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